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Abstract 

Random point fields which describe gases consisting of para-particles of any 
order are given by means of the canonical ensemble approach. The analysis for 
the cases of the para-fermion gases are discussed in full detail. 



1 Introduction 

Where do the statistics of random point fields come from? We examine what kind of 
random point fields follow from the para-statistics of particles. 

In the previous paper jHj, the boson and/or fermion point fields were derived by 
means of the canonical ensemble approach. That is, quantum mechanical thermal sys- 
tems of finite fixed number of bosons and/or fermions in the bounded box in M*^ were 
considered. By taking the thermodynamic limit to the position distribution of con- 
stituents, random point fields for boson and/or fermion gases of positive finite densities 
and temperatures on M'^ were obtained. There, the method was applied to construct the 
random point fields which describe gases consist of para-bosons (resp. para-fermions) 
of order 2. In the recent proceeding article JT], the argument for para-particle gases of 
order 3 is developed. 

In this paper, we pursue the project to the general case: we apply the method to 
statistical mechanics of gases which consist of para-particles of any order j9 G N. We 
will see that the random point fields obtained in this way are those of a = ±l/p given 
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in jOj. Our main result in this paper is 

Theorem The random point field for gas of para-fermions (resp. para-bosons of low 
density) of order p is equal in law to the convolution ofp independent copies of the usual 
fermion (resp. boson) point field for any p G N. 

We use the representation theory of the symmetric group (cf. e.g. [3 ISl EI)- Its 
basic facts are reviewed briefly in Sect. 2, along the line on which the quantum theory 
of para-particles are formulated. We state our main results in Sect. 3. Sections 4 and 
5 devoted to the full detail of the discussions on the thermodynamic limits for para- 
fermions and a few remarks on those for para-bosons, respectively. Some discussions are 
given in Sect. 6. 

2 Brief review on Representation of the symmetric 
group 

We say that (Ai, A2, ■ ■ ■ , A„) G N" is a Young frame of length n for the symmetric group 
Sn if 



E 



\j = N, Ai ^ A2 ^ ■ • ■ ^ A„ > 0. 

We associate the Young frame (Ai, A2, ■ ■ ■ , A„) with the diagram of Ai-boxes in the first 
row, A2-boxes in the second row,..., and A^-boxes in the n-th row. A Young tableau on 
a Young frame is a bijection from the numbers 1, 2, • • ■ , A^ to the N boxes of the frame. 

Let M.^ be the set of all the Young frames for Sn which have lengths less than 
or equal to p. For each frame in M^ , let us choose one tableau from those on the 
frame. The choices are arbitrary but fixed. 7"^ denotes the set of all tableaux chosen 
in this way. The row stabilizer of a tableau T is denote by 71{T) , i.e., the subgroup of 
Sn consists of those elements that keep all rows of T invariant, and C{T) the column 
stabilizer whose elements preserve all columns of T. 

Let us introduce the three elements 

o-iT) = ,,^ ,^, > a, b(T) = ,, ^ ,^, > sgn(o")cr 

and 

^(^) = M E E ^g^(^)^^ = ^Ta{T)b{T) 

of the group algebra C[5iv] for each T G 7^^, where dx is the dimension of the irreducible 
representation oiSN corresponding to T and ct = dx^ 4^7l(T)^C(T)/N\. As is known, 

a(Ti)a6(T2) = b{T2)aa{T,) = (2.1) 



hold for any a G Sn if ^2 — ^ ^i. The relations 

a{Tf = a{T), b{Tf = b{T), e{Tf = e{T), e{T,)e{T^) = Q {T, ^ T,) (2.2) 

also hold for T, Ti,T2 G T^^ . For later use, let us introduce 

d{T) = e{T)a{T) = CTa{T)b{T)a{T) (2.3) 

for T G T^ . They satisfy 

d{Tf = d{T), d{T,)d{T2) = (Ti^Ta) (2.4) 

as is shown readily from ()2.2j) and ()2.1|) . The inner product ( ■ , ■ ) of C[iS7v] is defined 
by 

(a, r) = (5o-r for a,T E Sjsi 

and the sesqui-linearity. 

The left representation L and the right representation R of Sn on C[iSAr] are defined 

by 

L{a)g = L{a) ^ g{T)T = ^ ^(r)(Tr = ^ ^(CT"V)r 

Tg5jv t£5]v t£Sn 

and 

/2(a)(7 = i?(a) 5^ (7(r)r = J^ ^(r)ra-i = J^ (7(ra)r, 

respectively. Here and hereafter we identify g : Sn — > C with X]re5 fi' (''")''" ^ C[5Ar]. 
They are extended to the representation of C[5Ar] on C[5Ar] as 

cr^T (7 T 

and 

RU)9 = 9f = Y.9(a)f{r)ar-' = $^ ( $^^(ar)/(r))a, 

where / = E. f{r)r = Er fi^'')^ = Er fi^)^''- 

The character of the irreducible representation of Sn corresponding to the tableau 
T G Tp^ is obtained by 

Xt{ct) = J2 {r,L{a)R{e{T))T) = J^ {r^crrSr)). 

We introduce a tentative notation as in ^ 

Xsi^) = E {r^Lia)Rig)r) = J^ {r,crrr')gil) = Yl dir-'ar) (2.5) 

for 5- = ErS-l^)^ e C[5Ar]. Then xt = Xe(T) holds. 



We consider representations of Sn on Hilbert spaces. Let Hi be a certain L^ space 
which will be specified in the next section and (d^Tii its A^-fold Hilbert space tensor 
product. Let U be the representation of Sn on ^^Hl defined by 

U{a)ipi (g) ■ ■ ■ o (^jv = V5a-i(i) ® ■ ■ ■ ® ^a-^N) for v^i, ■ ■ ■ , v^TV e Hl, 

or equivalently by 

(f/(o-)/)(xi,--- ,xjv) = /(x<^(i),--- ,a;^(7V)) for f e ^^Hl- 

Obviously, U is unitary: U{a)* = U{a~^) = f/(cr)~^. We extend U for C[iSAr] by linearity. 
Then U{a{T)) is an orthogonal projection because f/(a(T))* = U{a{T)) = U{a(T)) and 
dZl- So are f/(6(T))'s, t/(rf(T))'s and 

^PB = I] f/(rf(T)). (2.6) 

PpB is the projection operator to the subspace for para-bosons of order p. Note that 
Ran[/((i(T)) =Ran[/(e(T)) because d{T)e{T) = e(T) and e{T)d{T) = d{T). 

For para-fermions, we consider the transposed tableau T' of T G 7"^ by exchanging 
the rows and the columns of the Young tableau T. The transpose A' of frame A is defined 
in the same way. Then T' lives in A' if T lives in A. Clearly 

C{T') = 7^(T), 7^(T') = C{T) (2.7) 

and we define 

PpF = Yl f^(^(^')) (2.8) 

which is the projection operator to the subspace for para-fermions of order p. 

3 Para- statistics and Random Point Fields 

3.1 Para-fermions of Order p gN 

We first consider the quantum system of A^ para-fermions of order p in the box 
Al = [—L/2, L/2Y C M*^. We refer the literatures E El for quantum mechanics of 
para-particles. (See also |lj.) These literatures indicate that the state space of our 
system is given by Tif, ^ = Ppp ^^'Hl, where Tii = L^(Ai) with Lebesgue measure is 
the state space of one particle system in A^. We need the heat operator Gl = e^^'^ in 
Al, where A^ is the Laplacian in A^ with periodic boundary conditions at dA^. Then 

d 

spec Gl = {exp[— /^^J?!]; Qj G 27rZ/L}, 



i=i 



1 



4,TtGl = ^ 5^ exp (-/3(27r|n|/L)2) 



It is obvious that there is a CONS of Ti.^ jy which consists of the vectors of the form 

U{d{T'))ipl.J ® ■ ■ ■ ® (/'fc , which are the eigenfunctions of ®^ Gl- Then, we define a 
point field of N free para-fermions of order p in the box as in section 2 of '^ and its 
generating functional is given by 

^,, . _^,,^. ^ Erer/Tr^.^J(0^G)[/(d(rO)] 
^'^ L' J Erev T^ ^^«- [(®^G)f/(rf(r'))] ' 

where / is a nonnegative continuous function on Kl and Gl = G/ e^-^G/ . We first 
prove: 



Lemma 1 






'^TeT/ Ia^ detT'{GL{xi, Xj)}i<^ij<:Ndxi ■ ■ ■ dxN 



(3.2) 



J^TeT'^ Ja'^ detr/{Gj;^(a;j, a;j)}i^jj^7vrfa^i ■ ■ ■ dx^ 

Remark 1 : The state space H^ j^ = Ppp ®^ Ti-i is determined by the choice of 

the tableaux T's. The different choices of tableaux give different subspaces of ®^7iL- 

However, they are unitarily equivalent and the generating functional given above is not 

affected by the choice. In fact, Xt(c") depends only on the frame on which the tableau 

T is defined. 

Remark 2 : det^'^ = J2aes XT'{cr) Y[i=i ^ia(i) in ()H.2j) is called immanant. 

Proof : These expressions are derived by the relations 

Tr ^.n, [{®''GL)U{d{T'))] = Tr ^^n, [(®^Gz.)t/(e(T'))] 

and 

J2 X5MTr«iv^J(®^Gi)t/(a)] =iV!Tr^^.^J(®^Gi)f/(^)], (3.3) 

o-e5jv 

with g = e(T'). These relations can be shown by the use of ()2.5p . the cyclic property of 
the trace and the commutativity of U{t) with ®^G. For details, see 0. □ 

Now, let us consider the thermodynamic limit 

L, N ^oo, pl = N/L'^ ^ p > 0. (3.4) 

In the following, / is a nonnegative continuous function on M.'^ which has a compact 
support, and is fixed through the thermodynamic limit A^ /" W^. We identify the 
restriction of / to Al as / in Lemma ^ We get the limiting random point field on W^. 



Theorem 2 The finite random point fields for para-fermions of order p defined above 
converge weakly to the point field whose Laplace transform is given by 

EP^[g-(/,C>] = Det[l - x/l - e-/r,G(l + r.G)- Vl - e^f]" 
in the thermodynamic limit ^3-41 ), where r* G (0, oo) is determined by 
p f dp r,e-/^H (^rM^+nG)-'){x,x). 



p J {27:^1 + r^e-P\p\' 



3.2 Para-bosons of Order p gN 

We next consider the quantum system of A^ para-bosons of order p in the box A^, = 
[-L/2, L/2Y C R'^. It is obvious that there is a CONS of Hl^j^ which consists of the 

eigenfunctions for ®^Gl of the form U {d{T))Lp\,_ ■ ■ ■ (S> v^). • Then, we define a point 
field of N free para-bosons of order p as in section 2 of j|9j and its generating functional 
is given by 

ppB r -(/,0l ^ "^"^ ®^Hl[{®'^ Gl)Ppb\ 
^'"^^ J Tr^.^J(®^G^)P,^]' 

where / is a nonnegative continuous function on A^;, and Gl = G^ e^-'G^ . Then, we 
have: 

Lemma 3 

^,B r -urn ^ ^T,rr T..,s. XT(a)Tr,.^,[(^^G,)[/(a)] 



YliTi^r^ Ian detriGiixi, Xj)}dxi ■■■dx 



N 



Y^Ter^ Jj^N detT{GL{xi, Xj)}dxi ■ ■ -dx^' 



(3.6) 



p 



We again consider the thermodynamic limit ()3.4p . We get the limiting random point 
field on W^ for the low density region: 

Theorem 4 The finite random point fields for para-bosons of order p defined above 
converge weakly to the random point field whose Laplace transform is given by 



EPB[e-(/,0] = Det[l + ^/l^^r,G{l - r^G^Wl-e-f] " 
in the thermodynamic limit, where r* G (0, 1) is determined by 

p f dp r.e-^H^ (r,G(l-r.G)-i)(x,a;), 



p J (27r)'^ 1 - r,e-^H' 



p f dp e 

< Pc = 



-fi\p? 



P Jr^ (2vr)'^ 1 - e-/3bl= 



Remark : The high density region p ^ ppc is related to the Bose-Einstein condensation. 
We need a different analysis for the region. See [Tn| for the case of j9 = 1 and 2. 

4 Proof of Theorem |21 

It is enough to show the convergence of the generating functionals. In the rest of this 
paper, we use the results in jH] frequently. We refer them as, e.g.. Lemma 1.3.2 for Lemma 
3.2 of [9J. Although those in jH] are results for p = 1, their arguments hold for general 
p G N with obvious changes. Let ipT be the character of the induced representation 
Ind^^^Jl], where 1 is the one dimensional representation TZ{T) 3 a ^ 1, i.e., 

M^) = E {r,Lia)R{a{T))T) = Xa(T)(a). 

re5jv 



Since the characters xt and ipT depend only on the frame on which the tableau T lives, 
not on T itself, we use the notation x\ a^id t/'a ( A G M^ ) instead of xt and ipT, 
respectively. 

Let 6 be the frame (p— 1, ■ ■ ■ , 2, 1, 0) G M.^ . Generalize ipfj, to those /i = {pi, ■ ■ ■ , pp) G 
Z^ which satisfies Yl^=i f^j = N hj 

and 

•■p 1 



i'n = "ipwfi for p E Z^ and n E Sp such that np E M^ 



where Z+ = {0} U N. Then the determinantal form j2j can be written as 



Xx=Yl sgnn ijjx+s-ns- (4.1) 

TT€Sp 

Let us recall the relations 

Xr'(o-) = sgnaxr(o-), 'frier) = sgnaipria), 

where 

^T'M = $^(r,L(a)i?(6(T'))r) = xkH^) 

T 

denotes the character of the induced representation Ind^^,Jsgn], where sgn is the 
representation C{T') = TZ(T) 3 a ^-^ sgn a. Then we have a variant of ()4.H) : 

XX' = X^ sgnn ip X' +5' -{ttSY- (4.2) 

Now let us consider the denominator of p.ip . Let T G Tj^ live on /i = (/ii, • • • , pp) G 
Mp^. Thanks to (Q for g = b{T'),we have 

5^ ^T'(^)Tr^A^wJ(®^Gi)f/(a)] =iV!Tr^.^J(®^Gi)f/(6(T'))] 



where An = Ylres sgn(r)f/(r)/n! is the anti-symnietrization operator on &'1-Ll- In the 
last step, we have used 



KT) = n E 



sgno" 



-(T, 



where TZj is the symmetric group of yUj numbers which he on the j-th row of the tableau 
T. Now (Q yields 



^XA'MTr^.v^,[(®^Gz.)f/(a 



a"65]\ 



^ sgnvr ^ v'v+5/_(^5)/(a)Tr^iv^^[(®^Gi)[/((T)] 

AT! ^ sgn7rnTr^.,,-,,.(,)„J(®^^-^-+-(^-)Gi)A,^_,+.(,)]. 
7i"e5p i=i 



(4.3) 



Here we understand that Tr ^-it^l ((®"G'l)74„) = 1 if ra = and = if n < in the last 
expression. Applying the Cauchy integral formula to 

oo 

Det[l + ;2J] = J]z"Tr55n^[(®"J)A„] 

n=0 

where J is a trace class operator, we obtain that 



Tr^.^[(®"Gz.)A„] = / -^^Det[l + zGl], 



(4.4) 



where Sr{() = {z G C; |-2 — CI = '''}• Note that r > can be chosen arbitrary and 
the right hand side equals 1 for n = and for n < 0. Then we have the following 
expression of the denominator of ()3.ip : 

E E ^A'(a)Tr^.^J(®^G,)f/(a)] 



xeM^ o-e-S]v 



Det{l + ZjGL)dzj 

j+7T{j) + l ■ 



Ap{zu ■■■ , Zp) [n?=i Det(l + ZjGL)dzj] 



AT-E 



AeA4 



Sr{0)P 



uu 2^^^ 



(4.5) 



where Ap{zi, ■ ■ ■ , Zp) is the Vandermondian given by 



Ap{zi 



px'^L} } '^p; 



n 



yZi Zj^ 



^^i<j^P 



p—1 p— 1 

Z-l Zr, 



Z\ Z2 

1 1 



yP-\ 



(4.6) 



In the following, we simply write Ap{{z}) for A{zi, ■ ■ ■ ,Zp) when there is no danger of 
confusion. 

To make the thermodynamic limit procedure explicit, we take a sequence {LjvJArgN 
which satisfies N/Lf^- ^ p as N ^ oo. In the following, we set r = r^ E [0, oo) to be 
the unique solution of 



where 



Tr 



k 



1 + tGl, 



k 



N p-1 
,7^ 2 



(4.7) 



(4.8) 



is the averaged length of the rows in the Young tableau and [ ■ J represents the integer 
part. 

The existence and the uniqueness of the solution follow from the fact that the left- 
hand side of ()4.7|) is a continuous and monotone function of r. For details, see Lemma 
1.3.2 [for a = -1]. We also put 



Vn = Tr 



tnGi 



L(l + r^GL^)2 



(4.9) 



We will suppress the A^ dependence of vjy and L^- Since r^ —>■ r^ in the thermodynamic 
limit, we have k/{2 + r^) ^ v^ ^ k for large enough A^. See Lemma 1.3.5. [There r^r 
and r^ are written as zjq and z^, respectively.] 

Put 

,r^^ A^!Det[l + rG'jP 



^)p(v/^)i+p(p-i)/2r^ 



Then we have: 



Lemma 5 



lim Jp = —, 

N^oo p] 



I \A{y,,... ,yp)\5(J2y^\{e-'^"dy,>(). (4.11) 



3=1 i=i 



Proof : We set 
Then we have 



Uj = \j + p — j — k. 



p{p - 1) 



5Z^. = ^^0 = ^^ + 



pk e [o,p). 



i=i 



ui > z/2 > ■ ■ • > z/p ^ —k. 



The parametrization 



(4.12a) 
(4.12b) 
(4.12c) 



Zj = rexp[iXj/\'v 



(j = I,--- ,P) 



yields 



Det[l + ZjGl] = Det[l + rG^lDet [l + {zj - r)GL(l + tGl)-'] , 

Zj — r = r(e"'^/^ — 1) = r(^ism{xj/y/v) — 2sm'^{xj/2y/v)^, 

dzj = ire^^^^^dxjl ^^ 



and 



Jx) 



5^ n / ^ )A,({e^^/v^})(v^)i-^+^(^-i)/2 



p 



L / r- \ 111 / 






1 + (e* 



,/v^ 



-11 



rG, 



l + rG-L 



where the summation on v is taken over all v satisfying ()4.12|) . 
We consider two regions of x G (— a/wtt, a/wtt] 

1. small X region: \x\ ^ v^^^"^ ^ 

2. large a; region: |x| > v^l^"^ . 
In the large field region, we have 

Gl 



Det 



1+ z-r 



= Det 
^ Det 
^ exp ( — 



l + rGL_ 

X N tGl 



l-4sin^(-^) ^t. U-- 



rG/ 



+ rG 



L^ -I 



X X rG; 



sm , ^ , 

1 + r ^2/u^ 



1 + r 



sm 









(2^)T^rT^)<-p(— '^'"). 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 



(4.14) 



using ^ Gl ^ 1 and v = Oilz) = 0{N) and the boundedness oi r = r^ > uniformly 
in A^. 

In the small field region, we have the convergent expansion 



Det 



l + {z-r) 



G, 



l + rGri 



exp (e*^/v^-l)A;--(e"/^-l)^(A;-t;) 
+ ^LJ^(eWv^_l)^Tr 



/ Q lc"ir 1 

(l + 2_, ^f-^^ + Rn{x)) exp y—= - -x' 
e=3 ^^ 



rG, 



l + rG, 



(4.15) 



10 



where 



and 



IqI < const iV-^/*^ 

1 1 p 1 1 - 



3, ■ ■ ■ ,n — 1 



(4.16) 



sup |i?„(x)| =0(iV-"/^2^ 

|x|^l)l/12 



hold. We put 

y^cex = 5{x). 
e=3 
We choose n in ()4.15|) so large that 

Vf^^)^"p+p(p"^^/^llA II Hi? II =o(i) 



n-1 



(4.17) 



holds, i.e., n > 3{p — l)(j9 — 2). 

These arguments show that it is enough to consider the contribution from the small 
X region, and we have 



I' ^ j=i 



xA,i{e'- 



E n^+^^v/^^ 



l-p+p(p-l)/2 



d 



o{l), 



dui 



dx 



L ^-i^j^j/\^-^]/'^ 



^ 



u ^ j=l 

xAp({e^"/v^})(v^)i-P+P(P-i)/2| ^ ^(^^_ 
Thanks to the multi-linearity of the determinant Ap, we have 



(4,18) 



n 



/2^ 



iiyja;'j/v^-a;j/2 W C|g«^/v^|~) 



det < / ^=^e 
det|e-(^-'^^)'/2^ 



'^^i ^i{i~r.,)xj/v^-a;2/2 



p - 1 ^ f ^ 

1 ^ J ^ p 



p - 1 ^ l^ 

1 ^ i ^ P 

P ,,2 P-1 (/2 

A(K'"})-P(-E^-E^ 

i=l £=0 



(4.19) 



Since (0;)P(P-^)/2/^p({e^/''}) = Ap({v/i^(e''/'' - 1)}), the summation over all u satisfying 
the condition ()4.12|) yields 

lim y^(^)i-p+p(p-i)/2 X srm 



11 



= / 5iJ2yjWyi,---,yp)lle-'^/'dy,. (4.20) 

Jyi>--->yp j^i j^^ 

Here we have put i/j = Vjj \fv and regarded (\A^)^~^X]iy ^^ ^^^ integral of the suitable 
step function. Then, (a/w)^"^ ^^ — > / dyS(^j yj) is derived by the use of the dominated 
convergence theorem. The limit of the main term of ()4.18j) is given by ()4.2()|1 , which is 
equal to (J3J . We may see that the limit of the remainder vanishes from ()4.16p and the 
convergence of 



X 



ix/v^i w /-\l-p+p(p-l)/2 



Ap({e-/^n)(v^) 



<9"^ \ ,r^,l ..N { yr^V. 



V ,9 P-1 1,2 



= E(v^)"''(n(^)°'£i)^.({v^(-'""-i»)-p(-Ei-E 

— / '^(E%)(fll^)^^(?^i'---'^p)^"^^^^'^'^'(fl^^^)- (4-21) 

Jyi>->yp j=i j=i'^yj j=i 

We obtain this convergence by performing the differentiations in the second and the 
third members of ()4.21|) and applying the dominated convergence theorem. D 

The numerator is obtained just in the same way. That is, we replace Gl by Gl = 
G/ e~^Gl and introduce f = f^ and -u = t;Ar by 

rp rjyGiff J tnGlj^ 

Tr — = K, Tr ^r^^^ = V]^. 

l + r~NGL^ {l + r~NGL^y 

Gl — Gl is a positive operator of trace class such that Tr {G — G) = 0(||1 — e^-'^)!!!) 
and max {spec Gl} — max {spec Gl} = 0{L~'^). Then it follows from the definitions of 

rN,fN, vn and vj^ that 

^ f^ - r^ = 0{N-^), \v^ -vm\ = 0(1). (4.22) 

See Lemma 1.3.5 and Lemma 1.3.6 for detail. 

We define Jp similarly [v]^ is used instead of v^) and we get 

hm J, = ^,l W{yi.--- .y,)\5(Y,y^\{^-'^"dy,>Q 

by the very same argument as in the proof of Lemma El 



12 



Thus we have (writing r = r^, r = r^ and L = Ln 

\uex[i -|- tltlj/ 

Det 



Det[l +rGL]) Vf J UJ ^ 

Det[l+rGL^^'' 



Det[l+rGi] 



1 + ^^^^ZLg, 

l + rGr . 



r-f\^/i;\p(p-i)/4+i/2 J 
x(l + ^^) (-) -f- 



(i)' 



(4.23) 



Here 



Det[l + rGi] 



Det[l + rGL 
= Det 



Det 



1 + 



Det 



1- VI -e-/ 
1- vT 



1 + rGi 



(Gi - Gl) 



-f 



I + tGl 
r,G 



l-e~f 



VT 



=-/ 



■IP 



holds. For details, we refer Proposition 1.3.9 (and the argument on (c) in the proof of 
Theorem 1.3.1). The remaining factor of ()4.23p tends to 1 as iV — > cxd since v/v -^ 
1, Jp/Jp — * 1 and 



A^ log ( 1 + ^—^ ) - p log Det 



N^——^ -p 



r 



1 + i^^^G. 
l+rGi ■ 



r 



rGr 

Tr ^. + O(Ar-i) 



r 



^ 1 + rGi 

(A^ - pk) + 0{N-^) = 0{N-^). 



Finally, we also get 



P 



lim 



-Tr 



r^Gi 



p N^-^ L%" 1 + TnGl, 



{r,G{l + nG)-^){x,x) 



from (j4.7j) . (j4.8j) . N/Lj^ -^ p and Proposition 1.3.9. This completes the proof of Theorem 

121 □ 



5 Proof of Theorem SI 



rN 



In the case of para bosons, immanants (permanents) ^„(zSj^Y[j=i^L{xj,Xa{j)) and 

So-eSiv ni=i ^i(^j' ^cr{j)) appear in the denominator and the numerator, respectively. 
We can represent them by the Fredholm determinants by making use of Vere- Jones' 
formula USE]: 



f dz 

Tr ^.„[(®"G^)^„] = i -^_^Det[l - zG^ 

J SAO) ^^^^ 



-1 
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where Sn = ^res U{T)/n\ is the symmetrization operator on ®"'7iL and r G (0, 1) in 
this case. 

Using ()4.1|) , we get the following expression of the denominator of ()3.5|) : 



Y, $^XAMTr^.^J(®^Gz.)f/(a)] 



Ae>!^ fTG-Siv 






(5.1) 



where Ap{zi, ■ ■ ■ ,Zp) is the Vandermondian introduced in the previous section. 

We choose a sequence {LArJArgN which satisfies N/Lfj ^ p as A^ — >• oo. In this case, 
rjsf G (0, 1) denotes the unique solution of 

rGr 
Tr , ^^" =k (5.2) 

1 - rGLt, 



where 



as in (113) • We put 



,N p-1 , , , 

k=[j + '—l (5.3) 



Vn = Tr 



tnGi 



(5.4) 



l-rjvGL^)2. 

The remaining parts are almost the same as those in the para-fermion case. The reader 
may complete the proof of Theorem |3J following the previous arguments with the obvious 
changes. 

6 Discussion 

We have shown that 

the generating functional of the gas of para-fermions (resp. para-bosons of low density) 
of order p is equal to the p-th power of the generating functional of fermions (resp. 
bosons). 

The random point fields which we have obtained in this paper are a subset of those 
in [HI, where various properties of the point fields are examined. On the other hand, the 
authors of jB] obtained the point fields which do not follow from the representation theory 
of the symmetric groups which we discussed in this paper. Therefore it is interesting 
to consider physical interpretations of the point fields which do not follow from the 
representation theory of the symmetric groups. See e.g. P3]- 
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